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The  simpler  expression  for  the  median  noise  levels  (F^)  results  in 
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and  the  corresponding  values  computed  from  each  model.  The  most  accurate 
model  had  an  average  rms  error  for  all  times  and  seasons  of  1.25  dB;  whereas, 
the  least  accurate  had  an  rms  error  of  4.53  dB.  During  the  worst  season,  June, 
July,  August,  the  most  accurate  model  had  an  rms  error  of  1.4  dB;  whereas, 
tfye  least  accurate  had  an  rms  error  of  5.63  dB.  The  CCIR  noise  data  used 
to  generate  the  models  had  an  average  standard  deviation  of  the  error  at  1  MHz 
of  5.48  dB  and  an  error  of  6.67  dB  during  June,  July,  August;  the  numerical 
map  in  use  in  current  HF  prediction  programs  has  an  average  rms  error  of 
1.65  dB  and  an  error  of  2.02  dB  during  June,  July  and  August. 
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OBJECTIVE 

Develop  a  practical  but  simplified  atmospheric  noise  model  for  use  in 
estimating  system  signal-to-noise  ratios  in  micro-mini- computer  based  HP 
propagation  prediction  systems. 

RESULTS 

1 .  The  numerical  mapping  techniques  of  Zacharisen  and  Jones  were  used 
to  represent  worldwide  atmospheric  noise  at  1  MHz. 

2.  Each  four- hour  time  block  (in  universal  time,  UT)  for  each  of  four 
3-month  periods  was  represented  separately  rather  than  for  a  twenty- four  hour 
day,  resulting  in  24  numerical  maps. 

3.  This  numerical  representation  results  in  a  minimum  reduction  by  one- 
fifth  of  the  number  of  coefficients  in  computer  storage  at  one  time  over  that 
used  by  Zacharisen  and  Jones  for  their  model. 

4.  The  remaining  coefficients  would  be  stored  on  cassette  tape. 

5.  The  simpler  egression  for  the  median  noise  levels  results  in  sig¬ 
nificant  savings  in  computer  code  and  time. 

6.  The  models  produced  make  use  of  the  frequency  dependence  model  of 
earlier  representations  of  CCIR  Report  322. 

7.  The  average  rms  error  for  the  Zacharisen  and  Jones  model  for  all 
hours  and  seasons  was  1.65  dB;  the  corresponding  error  for  the  CCIR  Report  322 
data  at  1  MHz  was  5.4P  dB. 

8.  The  average  rms  error  for  all  times  and  seasons  for  a  model  version 
using  192  coefficients  per  time  block  was  1.25  dB. 

9.  The  average  rms  error  for  all  times  and  seasons  for  a  model  version 
with  110  coefficients  was  1.60  dB. 

10.  The  two  model  versions  with  77  and  55  coefficients  each  have  an 
average  rms  error  for  all  times  and  seasons  of  2.18  dB  and  2.94  dB, 
respectively. 
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11.  The  two  model  versions  with  50  and  35  coefficients  each  have  an 
average  rms  error  for  all  times  and  seasons  of  3.64  dB  and  4.53  dB, 
respectively. 

12.  The  least  accurate  season  for  all  model  versions  was  the  June,  July, 
and  August  season. 


RBCOMNKMDATIONS 

1.  Use  the  model  versions  with  192  coefficients  in  each  time  block  in 
either  Advanced  Prophet  applications  or  in  large  scale  HF  prediction  models. 

2.  Use  the  model  version  with  110  coefficients  in  each  time  block  in 
Advanced  Prophet  applications. 

3.  Use  the  model  versions  with  55  and  77  coefficients  each  in  each  time 
block  in  either  Advanced  Prophet  applications  or  Tektronix  4J52  mini-computer 
applications. 

4.  Use  the  two  model  versions  with  35  and  50  coefficients  each  in  each 
time  block  only  for  applications  for  which  the  memory  requirements  for  the 
mini-computer  demand  a  sacrifice  in  the  accuracy  of  the  model. 

5.  Use  Fourier  series  inte-roolation  to  obtain  Fm  at  1  MHz  for  an 
actual  hour  of  operation. 

6.  When  Fourier  series  interpolation  is  not  practical,  use  linear 
interpo 1 at ion . 

7.  At  HF  calculate  noise  values  for  galactic  and  man-made  noise  also. 

8.  Calculate  galactic  noise  only  for  frequencies  above  the  maximum 
expected  fQF2 . 

9.  Sum  the  noise  power  in  watts  of  the  three  noise  sources  and  convert 
back  to  dB  >1  watt. 

10.  Determine  the  reliability  of  the  combined  1  MHz  models  and  the 
frequency  dependence  model. 
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11 o  Determine  this  reliability  as  a  function  of  frequency,  time  of  day, 
geographic  location  as  well  as  for  the  maximum  longitude  and  latitude  harmonic 
nvmbers  of  the  1  MHz  numerical  maps. 

12.  Determine  the  accuracy  of  the  complete  Lucas  and  Harper  model  with 
the  same  data. 

13.  Use  measured  Fam  data,  measured  after  the  data  contributing  to  CCIR 
Report  322  were  measured,  as  well  as  data  used  to  develop  CCIR  Report  322. 
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is  for  any  type  of  business 
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Upper  decile  difference  from  the  median  noise  of  the  com¬ 
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Lower  decile  difference  from  the  median  noise  of  the  com¬ 
bined  noise 

Expected  value  of  the  signal  field- strength  required  for  a 
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Root-mean- square  noise  field  strength  for  a  1000  Hz  band¬ 
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logi0f) 

Effective  antenna  noise-factor  which  results  from  the 
external  noise  power  available  from  a  loss- free  antenna 
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Median  of  the  hourly  values  of  F  within  a  time  block 

cl 

Noise  factor  of  the  antenna  circuit  ( its  loss  in  available 
power) 

A  FORTRAN  program  used  to  determine  how  many  harmonics 
were  necessary  to  provide  an  adequate  fit  for  the  longi¬ 
tudinal  variation 

Noise  factor  of  the  receiver 

Noise  factor  of  the  transmission  line  ( its  loss  in  avail¬ 
able  power) 


galactic  noise  thermal  noise  that  originates  outside  the  earth  and  its 

atmosphere 
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(X) 


latitude  function  sin^X  cos^X 
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longitudinal  harmonic  maximum  (number  of  harmonics  N  = 
2J  +1) 
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Boltzmann's  constant  =  1.38X10-23 


joules  per  degree  Kelvin 
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Latitudinal  harmonic  maximum  (number  of  harmonics  N  = 
X  +  1  for  j  =  8,  9,  10;  for  0<j<7,  N  =  8) 
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Local  time 


man-made  noise  Industrial  noise 
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Combined  median  noise  from  all  sources 
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NBS  National  Bureau  of  Standards 

Ng  Median  galactic  noise 
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SI 


numerical  map  A  finite  series  of  mathematical  terms,  each  consisting  of 

a  product  of  a  numerical  coefficient  and  an  analytic  func¬ 
tion  of  the  geographic  coordinates 
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P  Received  signal  power  available  from  an  equivalent  loss- 

free  antenna  (dB) 
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A  numerical  coefficient  representing  atmospheric  noise  for 
the  hour  h  at  a  latitudinal  harmonic  k  and  a  longitudinal 
harmonic  j 

Signal-to-noise  power  ratio  required  (R  *  10  log1Qr) 
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residual  Predicted  noise  minus  observed  noise 
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Root  Mean  Square 
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or  similar  pursuits,  and  dwellings  are  no  more  than  one 
every  five  acres 

Ta  Effective  antenna  temperature  in  the  presence  of  external 
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TQ  Reference  temperature,  taken  as  288*  Kelvin 

t^  Hour  angle  in  degrees 

UT  Universal  time 


Voltage  deviation;  the  ratio  (dB)  of  the  rms  voltage  to 
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am 

0400,  0800,  1200,  1600,  2000  at  (X  ,0  ) 
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of  F^  using  the  Fourier  coefficients  generated  by  DATA 
FOUR 


A  FORTRAN  computer  program  that  calculates  the  coeffici¬ 
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dimensional  representation  produced  by  UT  NEQSOLj  it 
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IHTHOUUCTIOH 


The  determination  of  the  minimum  signal  level  required  for  satisfactory 
radio  reception  in  the  absence  of  other  unwanted  radio  signals  necessitates  a 
knowledge  of  the  noise  with  which  the  wanted  signal  must  compete.  There  are  a 
number  of  types  of  noise  that  may  influence  reception.  *fhese  can  be  divided 
into  two  main  categories  depending  on  whether  the  noise  originates  in  the 
receiving  system  or  externally  to  the  antenna.  External  noise  can  be  divided 
into  several  types.  The  most  usual  types  are  of  atmospheric,  galactic  and 
man-made  origin.  Below  about  30  MHz,  atmospheric  noise  usually 
predominates.  The  internationally  accepted  method  of  predicting  atmospheric 
noise  is  outlined  in  International  Radio  Consultative  Committee  (CCIR)  Report 
3221. 

The  purpose  of  this  report  is  to  describe  the  development  of  a  practical 
but  simplified  atmospheric  noise  model  for  use  in  estimating  system  signal-to- 
noise  ratios  in  micro-mini- computer-based  hf  propagation  prediction  systems. 
Previous  numerical  representations  of  atmospheric  noise  include  those  of  Lucas 
and  Harper  and  of  Zacharisen  and  Jones  .  The  Lucas  and  Harper  representation 
at  1  MHz  consists  of  24  numerical  maps:  one  for  each  of  six  4-hour  local  time 
blocks  and  for  each  of  four  3-month  periods.  (A  numerical  map  is  a  finite 
series  of  mathematical  terms,  each  consisting  of  a  product  of  a  numerical 
coefficient  and  an  analytic  function  of  the  geographic  coordinates.)  This 
model  uses  464  coefficients  for  each  time  block  and  2784  coefficients  for  a 
single  3-month  period.  The  frequency  dependence  is  represented  by  a  power 
series  of  least  squares  fit  to  the  frequency  dependence  curves  in  CCIR  Report 
322.  The  Zacharisen  and  Jones  model  is  a  worldwide  representation  of  1  MHz 
atmospheric  noise  in  universal  time.  The  entire  day  is  represented  with  one 
set  of  coefficients,  and  each  3-month  period  is  represented  with  only  960 
coefficients. 

The  approach  taken  to  develop  an  atmospheric  noise  model  suitable  for  use 
on  a  micro-mini-computer  was  to  try  to  reduce  the  size  of  the  numerical  map 
representing  atmospheric  noise  at  1  MHz  and  to  represent  each  four-hour  time 
block  separately  rather  than  for  a  twenty-four  hour  day.  Then,  assuming  that 
the  micro— computer  would  have  a  cassette  tape  unit,  the  coefficients  for  each 
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four-hour  time  block  for  each  season  would  reside  on  a  cassette.  Several 
versions  would  be  produced  -  each  with  a  different  number  of  coefficients  in 
each  time  block.  The  accuracy  of  these  versions  was  to  be  determined  as  a 
function  of  the  number  of  harmonics  in  the  representation.  A  user  would 
choose  the  version  with  an  error  he  was  willing  to  accept  consistent  with  the 
memory  space  available  in  his  micro-computer. 

The  numerical  mapping  techniques  of  Zacharisen  and  Jones  were  used  to 
represent  worldwide  atmospheric  noise  at  1  MHz.  The  data  used  to  model 
atmospheric  noise  was  the  same  as  used  to  prepare  the  local  time  contour  maps 
in  CCIR  Report  322 1  and  is  the  same  as  used  to  produce  the  two  previous 
representations  of  1  MHz  atmospheric  noise.  The  versions  produced  had  their 
worst  error  during  the  summer  months  of  June,  July  and  August.  Six 
atmospheric  noise  versions  were  produced:  (1)  a  numerical  map  with  192 
coefficients  having  an  average  rms  error  overall  of  1.25  dB  and  of  1.40  dB 
during  summer;  (2)  a  numerical  map  with  110  coefficients  having  an  average  rms 
error  overall  of  1.60  dB  and  of  1.76  dB  during  summer;  (3)  a  numerical  map 
with  77  coefficients  having  an  average  rms  error  overall  of  2.18  dB  and  of 
2.80  dB  during  summer;  (4)  a  numerical  map  with  55  coefficients  having  an 
average  rms  error  overall  of  2.94  dB  and  of  3.57  dB  during  summer;  (5)  a 
numerical  map  with  50  coefficients  having  an  average  rms  error  overall  of  3.64 
dB  and  of  5.25  dB  during  summer;  and  (6)  a  numerical  map  with  35  coefficients 
having  an  rms  average  error  overall  of  4.53  dB  and  of  5.63  dB  during  summer. 

The  material  presented  in  this  report  is  limited  to  the  development  of 
the  1  MHz  model  and  the  determination  of  its  accuracy.  Although  the  frequency 
dependence  model  is  included  in  an  appendix  for  completeness  the  treatment  of 
the  reliability  of  the  combined  1  MHz  model  and  the  frequency  dependence  model 
is  left  to  a  subsequent  report.  The  body  of  the  report  consists  of  five  main 
parts:  (1)  the  methods  of  atmospheric  noise  prediction;  (2)  a  description  of 
parameters  used;  (3)  the  numerical  mapping  process;  (4)  determination  of  the 
accuracy  of  the  models;  and  (5)  a  discussion  of  the  results.  Included  in  the 
appendices  are:  (1)  world  contour  maps  of  1  MHz  atmospheric  noise;  (2)  the 
frequency  dependence  model;  (3)  a  FORTRAN  computer  subroutine  of  the  1  MHz 
model;  and  (4)  a  computer  program  for  reading  and  printing  the  atmospheric 
noise  numerical  coefficients. 
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METHODS  OF  ATMOSPHERIC  NOISE  PREDICTION 


One  of  the  first  methods  of  worldwide  noise  prediction  was  that  published 
In  National  Bureau  of  Standards  (NBS)  Circular  No.  462. 4  in  this  publication 
noise  grade  maps  and  prediction  curves  were  given,  Indicating  the  noise  level 
in  terms  of  the  minimum  required  signal  strength  to  assure  radio-telephone 
communication  for  90%  of  the  time  in  the  presence  of  atmospherics. 

A  more  recent  publication,  NBS  Circular  No.  5575,  presented  the  same 
noise  grade  maps  as  used  in  the  previous  publication.  However,  the  prediction 
curves  had  been  revised  to  show  the  expected  median  levels  of  radio  noise 
during  each  season  instead  of  required  field  strength  for  90%  intel¬ 
ligibility.  The  method  of  interpreting  the  earlier  predictions  had  not  been 
entirely  clear,  and  the  new  method  of  presentation  was  used  to  remove 
ambiguities. 

An  even  more  recent  previous  prediction  method  was  contained  in  CCIR 
Report  65  (Revised),  which  was  adopted  by  the  CCIR  IXth  Plenary  Assembly,  Los 
Angeles,  1959. 6  The  information  presented  showed  the  expected  values  of  the 
average  noise  power  on  a  worldwide  basis  throughout  the  frequency  range  10  kHz 
to  100  MHz,  for  all  times  of  the  day  and  night,  and  all  seasons  of  the  year. 
The  presentation  was  revised  to  show  the  expected  median  levels  of  radio  noise 
during  four-hour  time  blocks  for  each  season.  The  parameter  presented  was  the 
median  hourly  value  of  the  average  noise  power  for  each  time  block,  and  the 
variations  in  this  parameter  showed  the  systematic  diurnal  and  seasonal 
variations  of  the  noise.  Variations  of  the  hourly  values  within  a  time  block 
were  treated  statistically,  and  their  extent  was  indicated  by  the  ratios  of 
the  upper  and  lower  decibel  values  to  the  median. 

Geographic  and  operating  frequency  distributions  of  atmospheric  noise 
from  CCIR  Report  No.  65  were  represented  numerically  by  Lucas  and  Harper  for 
use  in  computers. ^  Coefficients  were  tabulated  for  the  worldwide 
distributions  and  operating  frequency  dependence.  Representative  contour 
charts  of  all  distributions  were  shown  along  with  sample  longitudinal  and 
latitudinal  variations. 
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The  current  internationally  accepted  method  of  predicting  atmospheric 
noise  is  outlined  in  CCIR  Report  322.^  It  constituted  a  comprehensive 
revision  of  Report  65.  Although  the  original  form  of  presentation  was 
retained,  the  Report  was  improved  in  the  following  respects: 

-  the  original  predictions  were  modified  to  take  account  more  extensive 
and  reliable  data; 

statistical  information  was  presented  on  the  accuracy  of  the  modified 
predictions; 

the  fine  structure  of  the  noise  is  described  in  statistical  terms; 

-  each  time  block  has  its  own  frequency  dependence  curves  and  data  on 
noise  variability  and  character; 

-  methods  of  using  the  predictions  in  the  solution  of  operational  problems 
were  presented. 

Figures  1-3,  a  sample  of  the  figures  in  CCIR  Report  322,  show  the 
expected  median  values  of  atmospheric  noise,  in  dB  above  kTQb  for  summer, 

1200-1600  hrs  local  time.  The  figures  are  used  in  the  following  way.  The 
value  of  PM  for  1  MHz  is  found  directly  from  the  noise  charts  like  Figure  1 
for  the  time  block  (season  and  hour)  und***-  consideration.  Using  this  value  to 
select  the  noise  grade  curve,  the  value  of  Fam  for  the  required  frequency  is 
determined  from  the  curves  in  Figure  2.  The  parameters  op  ,  du,  and  are 
obtained  for  the  required  frequency  from  the  variability  curves  in  Figure  3. 
If  the  value  of  the  deviation  of  an  hourly  value  of  F  from  the  time  block 

Ot 

median  PM  (dB)  (D  »  Fa  -  PM) ,  or  the  value  of  op  is  required  for  any 
percentage  of  time  other  than  10%,  the  values  can  be  found  by  plotting  Du  and 
aDu  on  a  normal  probability  paper  (with  values  in  dB) ,  and  drawing  straight 
lines  through  O  dB  at  50%  and  the  10%  values.  Values  at  percentages  greater 
than  50%  can  be  determined  in  the  same  manner  using  and  oD^.  Plotted 
values  of  V^,  median  value  of  the  voltage  deviation  V^,  are  for  a  bandwidth 
of  200  kHz. 

Since  computer  methods  had  become  prominent  in  predicting  system 
performance  of  high  frequency  communication  circuits,  it  was  only  natural  that 
a  numerical  representation  of  CCIR  Report  322  would  be  made.  The  first 
representation  was  one  due  to  Lucas  and  Harper.2  Their  1  MHz  representation 
consisted  of  24  numerical  maps:  one  for  each  of  six  4-hour  time  blocks  (in 
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Figure  2 


Variation  of  radio  noise  with  frequency  (Suaner:  Figure  3 
1200-1600  h  LT) 


Data  on  noise  variability  and  character  (Su»er 
1200-1600  h  LT) 


local  time,  LT)  and  each  of  four  3-month  periods.  The  diurnal  variation  was 
obtained  by  linear  interpolation  between  the  numerical  map  values  taken  at  the 
central  hour  of  each  time  block.  Seasonal  variation  was  represented  by  a  step 
function.  Both  mean  and  frequency  dependence  variability  were  represented  by 
a  power  series.  Zacharisen  and  Jones  developed  a  representation  of  the  1  MHz 
atmospheric  noise  in  universal  time  (UT).3  They  reduced  by  two-thirds  the 
number  of  coefficients  necessary.  The  diurnal  variation  is  represented  by  a 
smooth  Fourier  series  (in  place  of  linear  interpolation),  and  the  seasonal 
variation  is  represented  by  linear  interpolation  (instead  of  a  step  function). 


DESCRIPTION  OF  PARAMETERS  OSED 


It  has  been  found  that  no  single  parameter  describing  the  noise 
environment  is  adequate  for  all  types  of  radio  service.  Thus,  a  number  of 
parameters  have  been  defined  to  evaluate  performance  of  various  systems.  This 
section  describes  noise  parameters  in  common  use  and  gives,  where  possible, 
the  relationships  between  them. 

The  most  useful  parameter  for  describing  the  noise  environment  is  the 
mean  noise  power  and  is  the  basis  of  noise  predictions.  The  noise  power 
received  from  sources  external  to  the  antenna  is  expressed  in  terms  of  an 
effective  antenna  noise  factor,  fa,  which  is  defined  by: 

fa  -  Pn/(kTob)  "  VTo  (1) 

where 

Pn  is  the  noise  power  available  from  an  equivalent  loss  free  antenna  in  watts; 
k  is  Boltzmann's  constant  (1.38x10-23  joules  per  degree  Kelvin);  TQ  is  the 
reference  temperature  288#K;  b  is  the  effective  receiver  noise  bandwidth  in 
Hz;  and  T#  is  the  effective  antenna  temperature  in  the  presence  of  external 
noise. 

Now,  using  capital  letters  corresponding  to  the  lower  case  letters  to 
denote  the  ratios  in  decibels,  Fa  represents  the  noise  power  density  per  Hz  in 
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decibels  above  kT  .  F  _  is  the  median  hourly  value  of  F  .  The  parameter  p„ 
o  am  a  n 

expressed  in  dB  above  1  watt  becomes: 

P  =  F  +  B  -  204  dBW  (2) 

n  a 

where  10  log1Q  kTQ  =  -204  and  B  =  10  log1Qb. 

The  assessment  of  the  performance  of  the  complete  receiving  system  can  be 
expressed  in  terms  of  its  operating  noise  factor,  f,  which  takes  into  account 
the  external  noise  as  well  as  noise  generated  within  the  receiving  system.  If 
it  is  assumed  that  the  receiver  is  free  from  spurious  responses  and  that  all 
elements  prior  to  the  receiver  are  at  the  reference  temperature  T  ,  then  f  is 
given  by: 


f 


fa  -1 


+ 

c  t  r 


(3) 


where  f  is  the  noise  factor  of  the  antenna  circuit  (its  loss  in  available 
c 

power) ;  ffc  is  the  noise  factor  of  the  transmission  line  ( its  loss  in  available 
power);  and  £  is  the  noise  factor  of  the  receiver. 

The  relation  between  the  noise  power,  pn,  available  from  a  loss  free 
antenna  and  the  corresponding  signal-to-noise  ratio,  r,  at  the  intermediate 
frequency  output  of  the  receiver,  can  be  derived  from  the  operating  noise 
factor  f  by: 


P 


n 


frkTQb 


(4) 


In  decibels  (dB),  P  expressed  in  dB  above  1  watt  becomes 
n 


P  =  R  +  F  +  B  -  204  dBW  (5) 

n 


where  F  is  the  noise  factor  f  in  decibels  above  kT0» 
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The  expected  value  of  received  power,  Pg,  required  for  a  particular  grade 

of  service  during  an  hour  when  the  antenna  noise  factor  is  F  ,  is  from 

8. 

equation  (5) 

Pg  =  +  R  +  B  -  204  dBW  (6) 

where  R  is  the  required  pre-detection  signal- to- noise  power  ratio  for  the 
given  bandwidth. 

In  some  cases,  it  is  desirable  to  deal  with  the  external  noise 
environment  in  terms  of  its  rms  electric  field  strength,  En.  For  a  short 
vertical  antenna  this  is  simply  related  to  the  noise  power  factor  by 

En  "  Fa  +  20  log1Q  fjjHz  -  65.5  (7) 

where  En  is  in  dB  above  1  uV/m  for  a  1  kHz  bandwidth  and  fMHz  is  frequency  in 
MHz.  The  value  of  the  field  strength  for  any  bandwidth  b  in  Hz  can  be  derived 
by  adding  (10  log-jg  b  -  30)  to  En-  Inserting  Fft  from  equation  (6)  in  equation 
(7)  gives  the  expected  value  of  the  field  strength  Eg  for  a  short  vertical 
rod: 

Ee  “  Pe  +  20  log10  fMHz  +  108,5  (dB  above  1  w/m>  (8) 

where  b  is  in  Hz. 
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NUMERICAL  MAPPING  PROCESS 


In  this  section  we  deal  with  the  problem  of  representing  the  worldwide 

geographic  variation  of  the  atmospheric  noise  parameter  F  at  1  MHz  for  a 

am 

fixed  hour  of  ijt  and  a  fixed  3-month  season.  The  procedures  described  here 
are  those  due  to  Zacharisen  and  Jones8.  These  procedures  were  based  on 
numerical  methods  for  representing  diurnal  and  geographic  variations  of 
Ionospheric  data  developed  by  Jones  and  Gallet.8 ' 1 0 

The  solution  to  the  problem  consists  of  well  defined  mathematical  opera¬ 
tions  which  have  been  programmed  for  use  on  large-scale  digital  computers. 
Input  to  the  computer  programs  are  values  of  FM  for  a  given  time  block  In 
local  time  and  a  3-month  period.  The  diurnal  variation  at  each  grid  point  Is 
represented  by  a  Fourier  analysis  of  the  local  time  block  values.  Then  a 

Fourier  analysis  is  made  of  the  longitudinal  variation  for  a  fixed  latitude. 
Next  the  worldwide  geographic  variation  of  each  Fourier  coefficient  Is 
expanded  in  a  series  of  functions  analogous  to  surface  spherical  harmonics. 
The  result  then  is  a  two  dimensional  representation  of  F^  In  latitude  and 
longitude  at  fixed  universal  times  for  each  season.  The  relatively  small 
table  of  coefficients  representing  F^  is  referred  to  as  a  "numerical  map." 


BASIC  INPUT  DATA 


The  values  of  F  ,  used  as  basic  Input  data  to  the  numerical  mapping 
programs,  were  the  same  as  those  used  to  prepare  the  LT  contour  maps  in  CCIR 
Report  322.  They  were  available  in  binary  form  on  magnetic  tape  for  each  of 
the  six  time  blocks  0000-0400,  0400-0800,  0800-1200,  1200-1600,  1600-2000,  and 
2000-2400  in  LT  and  for  each  of  the  3-month  seasonal  periods  December- 
February,  March-May,  June-August,  and  September-November . 

For  a  given  time  block  and  3-month  period,  the  values  of  F  were  given 

am 

at  8,400  grid  points  (X  ,  9  )  where 

m  n 


10 


100, 


(9) 


Xffl  =  90  -  1.8(m-1/2),  m  =  1,  2, 

.  360 

”  1/2),  n  —  1,  2,  •«•,  84*  (10) 

Here  X  denotes  latitude,  —90°  4X4  90° ,  and  9  denotes  longitude  measured 

eastward  from  0°,  0«4  6  4  3608.  The  grid  point  (X  ,  01)  is  at  89.1°N  and 

2. 1428571°E.  The  increment  in  latitude  is  1.8°,  and  in  longitude  it  is 

4.2857143°.  The  input  value  of  F__  at  (X  ,0  )  is  denoted  by  y 

am  m  n  ■*  Jm,n 

IOORIBR  RKPRKSENTATI ON  OF  D10HNAL  VAKIATIOH 


The  Zacharisen  and  Jones  numerical  representation  of  the  worldwide 
atmospheric  noise  was  in  universal  time  (UT) .  This  was  done  to  reduce 
inconsistencies  in  the  basic  input  data  near  the  poles  that  result  from  the 
use  of  local  time.  Greater  computing  efficiency  was  attained  with  the  UT  maps 
since,  for  most  applications,  calculations  are  made  for  fixed  instants  of 
UT.  To  develop  a  numerical  representation  in  UT,  it  was  necessary  to  convert 
the  data  from  local  time  to  universal  time. 

The  most  natural  method  for  representing  the  diurnal  variation  is  Fourier 
analysis,  since  atmospheric  noise  is  periodic  in  time.  Moreover,  the  trig¬ 
onometric  functions  associated  with  Fourier  analysis  are  automatically 
orthogonal  with  respect  to  equally- spaced  points  of  measurement;  hence,  the 
computational  problems  of  least  squares  fitting  are  greatly  simplified. 

The  basic  data  are  for  the  six  time  blocks  0000-0400,  0400-0800,  0800- 
1200,  1200-1600,  1600-2000,  2000-2400  in  local  time.  Zacharisen  and  Jones 
assumed  the  data  for  each  time-block  could  be  represented  by  the  central  hour, 
so  that  the  diurnal  variation  was  defined  at  the  6  hours  0200,  0600,  1000, 
1400,  1800,  2200,  in  LT.  We  let  y1 ,  y2,  ...,  y6  denote  the  six  time  block 
values  ya/J1  at  Xffl,  Qn  corresponding  to  the  hour  angles 

tt  -  15i  -  30,  i  -  2,6,.. .,22  (11) 

respectively.  By  Fourier  analysis  of  the  data  for  these  six  hours,  a 
continuous  diurnal  representation  is  given  by 
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Y(t  )  =  a  +  t  la  cos  jt  +  b  stnjt  }  +  a  cos  3t 
i  o  j=1  j  i  j  i  3  t 

where  the  Fourier  coefficients  are  given  by 


(12) 


ao  6  £  yh 

aj  =  7  6  yh  cos  jti 

bj  =  ?  £  yh  sin  jti 


j  j  =  1>2 


(13) 


Here  £  denotes  a  sum  over  h  =  1,  2,  ...»  6.  In  equation  (13) 

h 

that  the  initial  time  value  was  at  0000  rather  than  0200  hours, 
a  phase  shift  of  30  degrees  in  equation  (11). 


it  was  assumed 
This  required 


Figures  4-6  are  examples  of  such  a  Fourier  representation.  The  curves 
labeled  2.5  harmonics  were  obtained  using  equations  (11)  through  (13).  The 
curves  labeled  2.0  harmonics  are  for  the  case  when  a^  was  set  to  zero  or  only 
five  coefficients  were  used.  As  can  be  seen  six  data  values  (shown  by  O)  are 
represented  exactly  by  the  Fourier  representation  using  the  complete  set  of 
six  coefficients;  whereas,  the  set  of  five  coefficients  produce  values  as  much 
as  5  dB  off  at  the  original  data  points. 


By  means  of  the  Fourier  representation,  a  computer  program,  called  DATA 
FOUR,  determined  these  six  Fourier  coefficients  for  each  of  the  OT  hours  h  = 
0000,  0400,  0800,  1200,  1600,  2000  at  each  of  the  8400  grid  points  (*m,  ©n) 
for  each  season.  A  second  program,  called  UT  ANOIS,  computed  values  of  y 
F  for  each  of  the  UT  hours  h  -  0000,  0400,  0800,  1200,  1600,  2000,  at  each 
of  the  8400  grid  points  (^m,0n)  and  stored  them  on  magnetic  tape.  The  UT  grid 
points  are  coincident  with  that  for  the  LT  maps,  and  each  UT  grid  contains  84 
longitude  by  100  latitude  points. 
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(dB  ABOVE  kT  b|  F  (dB  ABOVE  kT  bl  F  „  (dB  ABOVE  kT  bl 

w  am  o  am  o 


Figure  4.  Diurnal  variation  of  basic  atmospheric  noise  data, 
June,  July,  August  (20.7°N,225.0°E) 


Figure  5 .  Diurnal  variation  of  basic  atmospheric  noise  data, 
June,  July,  August  (26.1°  S,  336.43°E) 


Figure  6.  Diurnal  variation  of  basic  atmospheric  noise  data, 
June,  July,  August  (0.9°N,272.14°E) 
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LONGITUDINAL  VARIATION  FOR  FIXED  LATITUDE 


The  first  step  in  representing  the  worldwide  geographic  variation  of  the 

atmospheric  noise  parameter  F  for  a  fixed  hour  of  UT  and  a  fixed  3-month 

season  is  to  hold  the  latitude  fixed  and  consider  the  representation  of  the 

periodic  function  of  longitude  0  defined  at  84  equally  spaced  points  0  . 

Zacharisen  and  Jones  Fourier  analyzed  the  data  y^h*  at  the  data  points  (1  , 

in*  n  111 

©n)  using  an  express^  n  of  the  form 


(h)  J  ,  (h)  „  (h)  , 

Z  (0)  =  a  +  I  [a  cosj  0  +  b  sin  j  ©],  (14) 

h,m  m,o  j=1  m,j  m,j 


where  J  denotes  the  number  of  harmonics  (2J  +  1  <  N  =  84). The  coefficients  are 
defined  by  the  relationships 


(h)  1  N  (h) 

a  =  —  £  y 

m,o  N  n=1  m,n 


(h) 

L 

m,  j 


2  N  (h) 

—  Z  y  cos  j  0  , 

N  n=1  m,n  n 


(15) 


(h)  2  N  (h) 

b  a  —  Z  y  sin  j  0  . 

m,j  N  n=1  m,n  n 


At  the  outset,  it  was  not  known  precisely  how  many  harmonics  would  be 
needed  to  provide  an  adequate  fit  for  the  longitudinal  variation,  zacharisen 
wrote  a  computer  program  called  FOUR  RESD  to  analyze  both  the  mapping  residual 
and  the  coefficients  themselves  to  aid  in  the  determination  of  the  number  of 
harmonics  to  be  used.  The  program  determines  the  standard  deviation  of  the 


mapping  residual  for  J  =  5  to  25  where  the  residual  equals  the  evaluation 

of  y^  at  each  longitude  using  the  Fourier  coefficients  minus  the  actual 
m#n  (h) 

value  of  y  .  The  quantity 
m,n 


e 

®,J 


2 


D 


[ 


m,  J 
84-2 J-1 


1/2 


(16) 


is  used  as  an  estimate  of  the  standard  deviation  of  the  residuals,  84-2 J-1 
being  the  number  of  degrees  of  freedom  after  subtracting  one  for  each  term  in 
equation  (14).  The  square  of  the  residuals  is  given  by 


2 

D 

m,  J 


84  ,  (h) 

1  « 

n“  1  m ,  n 


,2 

Z  (0  )] 

h,m  n 


(17) 


In  addition  it  calculates  the  average  Fourier  spectra  (means  of  100  values) . 
Fourier  spectra  are  calculated  using 


r  (hi  2  .  (h) .2  ,  (h) .2 

[Cl  -  [a  ]  +  [b  ]  . 

m,j  ® » j  ® ,  j 


( 18) 


The  results  indicated  to  Zacharisen  and  Jones  that  the  series  could  be 
terminated  at  J*10  resulting  in  21  coefficients  for  each  latitude.  Rarely  did 
they  feel  that  J  need  be  larger  than  17 . 

Having  determined  the  necessary  maximum  number  of  harmonics  J  for  the 

longitudinal  variation,  Zacharisen  wrote  a  computer  program,  called  UT  FOUR, 

to  calculate  the  coefficients  and  b^  in  equation  (14)  at  each 

®,  j  m,j 

latitude  Ajl  for  each  hour  h  and  each  season.  Fourier  coefficients  for 
harmonics  up  to  J>18  were  output  on  a  data  tape.  As  these  coefficients  are 
orthogonal  (the  elimination  of  one  does  not  affect  the  value  of  another),  the 
number  used  in  practice  can  be  reduced. 


15 


LATITUDINAL  VARIATION  OF  FOURIER  COEFFICIENTS 


Having  determined  the  longitudinal  variation  in  the  equation  (14) 
using  a^J  and  b£^  for  fixed  latitudes  X^,  the  latitudinal  variation  of  the 
Fourier  coefficients  themselves  can  be  represented  by  series  of  the  form 


L  (X) 

K 

(h) 

=  E 

P 

k=0 

K. 

(h) 

>  (X) 

=  z 

q 

h'j 

k=0 

G  (X),  j 


Of  1/  •  *  •  f  J8*  1 0 


G  (X)  , 


j  5  1  #  2 ,  •  •  •  t  *1  38  10  , 


(19) 


where  the  latitude  functions  G.  .  (X)  are  defined  by 

3 » K 


Gj,k  (X)  =  slnk  X  cosj  X  (20) 

Here  X  denotes  the  geographic  latitude  -90®<  X  <90*.  The  function  G.  (X) 

3  rk 

j  /  0  has  a  built-in  constraint  at  the  geographic  poles  (X  =  ±90®)  in  that  the 
cosine  factor  is  equal  to  zero.  This  property  resembles  that  of  the  classical 
Legendre  polynomials.  However  orthogonal  series  of  Legendre  tend  to  blow  up 
near  the  end  points.11  The  function  given  in  (20)  is  as  stable  at  the  end 
points  as  in  the  center. 

The  coefficients  an<j  aih!  were  determined  by  the  method  of  least 

squares  using  the  FORTRAN  program  UT  NEQSOL. J  The  estimates  of  standard 
deviations  of  residuals  Ej,*  and  Fj,k  9iv«n  aa  follows 
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100  r  (h) 
I  [a 

n^l  m, j 


a  (X  )-  /( 100-K-1 ) 

h,j  m 


1/2 


100  (  (h)  2 

l  [b  -  b  (X  )*•  / ( 1 00-K-1 ) 

m>*1  m,j  h,j  m 


1/2 


(21) 


were  used  to  measure  the  goodness  of  the  approximation.  Here  100-K-1  is  the 
number  of  degrees  of  freedom  remaining  after  subtraction  of  one  degree  for 
each  term  in  the  series  for  a  ^  j(X)  or  b^  ^(X).  Zacharisen  and  Jones  plotted 
Ej,lc  and  Pj,k  gainst  K  for  several  hours  and  seasons.  They  determined  that 
K~9  was  a  good  average  cutoff,  which  was  used  fox  harmonics  0  <  j  <  7.  For 
harmonics  j  =  8,  9,  and  10,  they  choose  the  cutoff  K  -  6.  This  results  in  192 
coefficients  per  hour  season  to  represent  8400  data  points. 

These  analyses  also  showed  that  an  average  cutoff  of  X>6  for  all 
harmonics  would  result  in  only  a  small  Increase  in  the  variance  ox  the 
residuals  and  would  provide  large  reductions  in  storage  of  coefficients  and  in 
computer  time  required  to  evaluate  the  numerical  map.  Using  K=6  results  in 
147  coefficients  per  hour  per  season  to  represent  F^.  However,  because  the 
functions  used  in  (19)  are  not  orthogonal,  the  whole  series  must  be  recomputed 
if  terms  at  the  end  of  the  series  are  to  be  cut  off. 


TWO-DIMENSIONAL  GEOGRAPHIC  REPRESENTATION 


Thus  far  we  have  described  a  representation  of  the  two-dimensional  geo¬ 
graphic  representation  of  Fam  at  1  MHz  for  a  fixed  UT  hour  h  and  season.  This 
representation  has  the  form 
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X  (X,  0) 
h 


G  (X) 
o,k 


K 

I 

k=0 


<h) 

P 

o,k 


J  f  K  (h) 

+  £  {  (  £  P  G  (X)  )  cos  j© 

j“1  k=0  j,k  j,k 


(22) 


K  (h) 

+  (  £  q  G  (X))  sin  j© 
k=0  j,k  j,k 


where  K=9  for  0<  j  <  7,  K=6  for  j  -  8,  9,  10,  and  J  =  10.  G ^  k  (X)  is  given 
by  equation  (20).  The  coefficients  Pjh£  and  ^j^k  are  tound  using  UT 
NEQSOL  as  described  in  the  previous  section  for  6  hours  of  UT,  h**0000,  0400, 
0800,  1200,  1600,  2000  for  each  3-month  season. 

Zacharisen  and  Jones  carried  the  mapping  progress  one  step  further.  The 
numerical  coefficients  from  UT  NEQSOL  were  Fourier  analyzed  in  UT  by  a  program 
called  UT  FOUR3,  producing  one  map  for  each  three-month  season.  They 
originally  used  six  coefficients  for  this  purpose,  but  later  dropped  the 
highest  order  time  harmonic  coefficient.  Their  final  map  was  192  x  5  *  960 
coefficients. 
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DETERMINATION  OF  MODEL  ACCURACY 


The  approach  taken  in  the  development  of  the  model  was  that  several 
versions  would  be  produced  -  each  with  a  different  number  of  coefficients 
representing  in  each  time  block.  Then  the  accuracy  of  these  versions  was 

to  be  determined  as  a  function  of  the  number  of  harmonics  in  the  representa¬ 
tion  (number  of  coefficients).  The  user  would  choose  the  version  with  an 
etror  he  would  be  willing  to  accept  consistent  with  the  memory  space  available 
in  his  micro-mini-computer , 

There  are  two  methods  for  assessing  the  accuracies  of  the  models.  The 

first  would  be  to  study  the  residuals  between  values  Ym^n  ° Fam  an<* 

corresponding  values  computed  from  the  numerical  map  using  equation  (22).  The 

second  would  be  to  produce  UT  contour  maps  of  F  in  dB  above  kT  b. 

am  o 

Examination  of  the  geographical  detail  for  each  of  the  versions  would  indicate 
differences  between  the  various  choices. 

In  the  comparisons  to  follow,  the  accuracies  are  given  in  terms  of  the 
longitudinal  harmonic  number  J  and  the  latitudinal  harmonic  number  K.  Table  1 
gives  the  number  of  coefficients  as  a  function  of  these  maximum  harmonic 
numbers. 


J 

NUMBER  OF  COEFFICIENTS 

K  =  9  K  =  6 

K  = 

10 

192 

147 

105 

9 

178 

133 

95 

8 

164 

119 

85 

7 

150 

105 

75 

6 

130 

91 

65 

5 

110 

77 

55 

4 

90 

63 

45 

3 

70 

49 

35 

2 

50 

35 

25 

1 

30 

21 

15 

Table  1.  Number  of  Coefficients  as  a  Function 
of  the  Longitude  and  Latitude  Maximum 
Harmonic  Numbers 
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MAPPING  RESIDUALS 


To  develop  a  sufficient  knowledge  about  the  capabilities  of  the  various 
versions  of  the  representation  of  F  possible,  residuals  were  computed 
between  values  y^*1^  of  Faffl  and  corresponding  values  computed  from  numerical 
maps  using  equation  (22).  The  term  residual  is  used  with  the  following 
meaning: 

residual  =  (predicted  value)  -  (observed  datum).  (23) 


Certain  statistical  measures  of  this  term  have  proved  useful  in  past 
studies  in  comparing  predicted  and  observed  data.  These  include  the  average 
residual  (av.  res.)  and  the  root  mean  square  residual  (rms  res.).  The  average 
residual  locates  the  center  of  the  distributions  of  error  and  is  often 
referred  to  as  the  bias  in  the  estimate.  The  root  mean  square  residual  is  a 
measure  of  the  dispersion  in  the  error.  In  fact,  the  rms  res.  is  the  standard 
deviation  of  the  error  about  the  origin  (zero  bias)  and  is  related  to  the 
standard  deviation  about  the  mean  according  to 


2 

a  =  v 

2 


2 

v 

1 


(24) 


where  is  the  mean  square  error  (the  square  of  the  rms  error)  and  is  the 
bias.  When  the  bias  is  small  or  nearly  zero,  then  the  standard  deviation  and 
the  rms  error  are  nearly  the  same.  Otherwise,  the  rms  error  is  larger  than 
the  standard  deviation. 

The  sampling  of  v(h)  consisted  of  values  at  a  rectangular  array  defined 

n 

by  the  intersections  of  42  equally  spaced  longitudes  and  30  equally  spaced 
latitudes  between  X  =  81  ®N  and  X  =  81° S,  or  about  one-sixth  of  the  original 
grid  of  100  latitude  by  84  longitude  data  points.  Since  residuals  at 
the  high  latitudes  |A|  >  81®  were  considerably  larger,  they  were  omitted  from 
these  computations. 


20 


For  comparison.  Table  2  gives  the  rms  error  for  Zacharisen  and  Jones 
model.3  Examination  of  the  table  reveals  that  the  rms  error  is  largest  during 
the  summer  months  of  June,  July  and  August.  The  largest  dispersion  in  the  rms 
error  with  time  occurs  during  the  spring  months  of  March,  April  and  May.  The 
bias  for  the  2iacharisen  and  Jones  model  is  for  all  practical  purposes  zero. 


Hour 

(UT) 

December 

January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

0000 

1.66 

1.35 

2.07 

1.43 

0400 

1 .53 

1.34 

2.11 

1 .45 

0800 

1.71 

1.57 

2.08 

1.57 

1200 

1.84 

1.74 

1.94 

1.69 

1600 

1.54 

1.42 

1.93 

1.49 

2000 

1.41 

1.30 

1.97 

1.28 

Average 

1.62 

1.45  . 

1.01 

1.49 

Std.  Dev. 

0.15 

0.17 

0.08 

0.14 

Table  2.  RMS  Residual  in  dB  of  the  Zacharisen 
and  Jones  Model. 


The  initial  comparison  was  for  the  maximum  number  of  harmonics  K  for  the 
latitudinal  variation  equal  9.  The  maximum  number  of  harmonics  J  for  the 
longitudinal  variation  was  varied  from  1  to  10.  The  bias  was  found  to  be 
practically  zero.  Figure  7  shows  how  the  average  rms  error  in  dB  varies  with 
longitudinal  harmonic  number  J;  the  average  is  the  average  over  the  six  hours 
of  the  day.  In  the  figure  winter  means  December,  January,  February;  spring 
means  March,  April,  May;  Bummer  means  June,  July,  August;  and  fall  means 
September,  October,  November.  For  J  equal  5  through  10  there  is  little 
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variation  In  the  rms  error  except  for  the  seasonal  difference.  Below  J=5 
there  is  a  rapid  increase  of  rms  error  with  decreasing  J.  In  fact  at  J=1 ,  the 
summer  has  an  rms  error  of  8.2  dB.  Figure  8  gives  the  standard  deviations 
with  time  of  the  rms  error.  The  standard  deviation  does  not  vary  appreciably 
for  J=5  to  J=10  but  rapidly  increases  with  decreasing  J  below  5. 

The  next  step  was  to  assess  the  effect  of  allowing  the  maximum  number  of 
harmonics  K  for  the  latitudinal  variation  to  vary.  Since  the  functions  used 
in  the  latitudinal  variation  are  not  orthogonal,  it  was  necessary  to  run  the 
computer  program  UT  NEQSOL  for  each  value  of  K  being  considered.  In  addition 
coefficients  previously  determined  for  K=9,  coefficients  were  also  calculated 
for  K=6  and  K=4.  These  were  then  used  to  determine  the  residuals  as  a 
function  of  longitudinal  harmonic  number  J  and  season.  In  all  cases  the  bias 
was  practically  zero.  Figures  9-12  present  the  average  rms  error  over  time  as 
a  function  of  longitudinal  harmonic  number.  Figure  9  is  for  the  winter  season 
of  December,  January,  and  February;  figure  10  is  for  the  spring  season  of 
March,  April  and  May;  figure  11  is  for  the  summer  season  of  June,  July,  and 
August;  and  figure  12  is  for  the  fall  season  of  September,  October,  and 
November.  In  all  of  the  figures,  there  are  two  apparent  effects.  The  first 
is  in  the  region  for  J=5  to  J=10.  In  this  range  the  variation  is  mainly  due 
to  changing  K.  The  effect  of  changing  J  is  not  as  pronounced  here  as  below 
J=5 .  The  change  caused  by  varying  K  is  not  consistent  from  season  to  season; 
the  change  is  most  pronounced  in  the  winter  and  summer  season.  Here  again  it 
is  not  linear.  During  winter  it  changes  slower  than  during  the  summer 
months.  Below  J=5,  even  though  the  average  rms  error  increases  with 
decreasi  »g  value  of  K,  the  main  effects  are  caused  by  decreasing  J.  This  is 
particularly  true  at  J=2  and  below. 
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Average  rms  error  In  dB  of  1  MHz  atmospheric  noise 
model  as  a  function  of  latitudinal  and  longitudinal 
harmonic  numbers,  September,  October,  November 


WORLD  NAP  CONTOURS 


Using  a  computer  evaluation  of  the  numerical  maps,  one  can  produce  UT 
contour  maps  of  F^  in  dB  above  KTQb  directly  on  a  plotter.  This  is  parti¬ 
cularly  useful  in  malting  choices  between  the  various  models.  Examination  of 
the  geographical  detail  will  reveal  missing  contours  of  importance.  The 
absence  of  contours  of  high  level  would  be  particularly  significant. 

Contour  maps  were  plotted  for  the  hours  0000  and  1200  during  the  seasons 
December,  January,  February  and  June,  July,  August.  A  contour  map  was  drawn 
for  each  harmonic  number  combination  of  J=1Q  and  K*9i  of  J=5  and  K”9,6,4;  and 
of  J=>2  and  K=-9,7.  A  contour  plot  was  also  drawn  for  each  hour  and  season  for 
the  2acharlsen  and  Jones  model.  Figures  13-19  are  for  hour”00Q0  during  June, 
July,  and  August.  Because  the  residual  analysis  showed  this  hour  and  season 
to  have  a  particularly  high  rms  error,  it  was  chosen  to  illustrate  the  effects 
of  reducing  the  number  of  coefficients.  A  contour  map  is  also  Included  for 
the  Zacharlsen  and  Jones  model  for  the  same  hour  and  season.  The  remaining 
contour  maps  are  presented  in  Appendix  A. 

Figure  13  is  for  the  Zachatisen  and  Jones  model,  and  Figure  14  is  for  the 
model  with  J=10,  K=9.  There  is  slightly  more  detail  in  Figure  14  in  the  equa¬ 
torial  region.  In  particular  there  is  a  contour  in  Figure  14  of  10  dB/kTb 
missing  in  Figure  13.  Also,  in  the  left  portion  of  Figure  13  for  the 
Zacharlsen  and  Jones  model,  the  detail  is  smoothed  out  more.  In  the  northern 
polar  regions  there  is  a  break  in  the  20  dB/kTb  contour  level  in  Figure  14 
missing  in  Figure  13.  Figure  15  for  J“5,  K“9  is  very  much  like  that  for  the 
Zacharlsen  and  Jones  model  except  it  has  a  little  more  detail  in  the  polar 
region. 

As  the  number  of  harmonics  is  reduced  further,  the  high  level  contours 
begin  to  disappear.  In  Figure  16  for  J*5,  K*6  two  high  level  contours 
disappeared  --  a  85  dB/kTQb  and  a  90  dB/kTQb  in  the  equatorial  region.  Also  a 
20  dB/kT0b  contour  level  curve  disappeared  in  the  northern  polar  region.  In 
Figure  17  for  J*5,  K»4  >  there  is  a  considerable  loss  of  detail  on  the 
righthand  in  the  equatorial  region.  Both  Figures  18  for  J-2,  K-9  and  19  for 
J«2  and  K-6 ,  show  a  loss  of  detail  over  the  whole  contour  plot  from  the 
previous  versions.  This  is  more  noticeably  true  in  the  latter. 
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420.0 


Figure  13.  «or:.d  contour  map  of  ranl  at  1  MH*  for  ZacharUan  and 
Jones  model  foi  Jane,  July,  August  0000  hrs  UT 


Figure  14.  World  contour  nap  of  rm  «t  1  MHr  for  a  aodel  with 

192  coefficients  per  tine  block  (J*10,  K*9)  for  June 
July/  August  0000  hrs  UT 


LUTITUOC  -»  T>  «M 


Figure  16*  World  contour  Mp  of  F  at  1  MHe  for  a  nodal  with  77 
coefficient*  per  time  block  {J-5,  K»6)  for  June, 


.Tnltf .  0000  hr*  ITT 


Figure  17.  Nor  Id  contour  tup  of  ?u  *t  1  mi  for  a  aodel  with  55 
coefficient*  per  tiae  block  (J-5,  K-4)  for  June, 

July  ■  .\ugu*t  0000  hr*  OT 


60.0 


420.0 


FlgUre  18.  World  contour  nap  of  Vmm  et  1  HHr  for  a  aodel  with  50 
coefficient*  per  time  block  (J-2,  K-9)  for  June, 

July,  kuguat  0000  hra  OT 
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DISCUSSION  OF  RESULTS 


This  report  presented  the  development  of  several  minicomputer  atmospheric 
noise  models.  This  was  done  using  the  Zacharisen  and  Jones  numerical  mapping 
techniques.  Several  versions  were  produced  —  each  with  a  different  number  of 
coefficients  representing  Fam  in  each  time  block.  The  accuracies  of  these 
versions  were  presented  in  the  previous  section  as  a  function  of  the  number  of 
harmonics  in  the  representation  (number  of  coefficients) .  The  user  then  can 
choose  the  version  with  an  error  he  could  accept  consistent  with  the  memory 
space  available  in  his  micro-mini-computer. 

It  is  the  purpose  of  this  section  to  provide  some  guidance  in  the 

selection  of  a  model,  to  review  methods  available  for  interpolation  between 

time  blocks,  and  to  describe  how  these  model  values  of  F  can  be  used  with 
9  axn 

values  of  Fam  produced  by  models  of  man-made  noise  and  galactic  noise. 

CHOICE  OF  MODEL 

In  addition  to  the  memory  space  available  in  a  particular  minicomputer, 
several  other  factors  would  influence  the  choice  of  a  noise  model.  These 
include  the  accuracy  of  the  CCIR  322  noise  data  used  to  generate  the  models, 
the  accuracy  of  the  Zacharisen  and  Jones  model,  and  the  accuracy  of  the  noise 
models. 

accuracy  of  CCIR  322  Holae  Data 

In  the  previous  section  the  accuracy  of  the  noise  model  was  determined  by 
determining  the  residuals  with  the  data  used  to  generate  the  models.  These 
data  were  the  same  as  those  used  to  generate  the  contour  plots  in  CCIR  Report 
322.  Hence  the  accuracy  of  theBe  data  has  some  bearing  on  the  choice  of 
models.  If  the  rms  residual  is  low  compared  to  the  error  of  the  CCIR  322 
data,  then  it  is  probably  an  acceptable  model.  However,  if  it  is  nearly  equal 
to  or  larger  than  that  for  the  CCIR  322  noise  data,  then  combined  errors  would 
be  a  better  Indicator  of  the  actual  error  in  the  model. 

The  error  in  the  values  of  F^  presented  in  CCIR  Report  322  is  given  for 

each  time  block  for  each  season  as  op  .  It  appears  as  a  curve  as  a  function 

am 

of  frequency  like  Figure  3. 
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December 

March 

June 

September 

Hour 

January 

April 

July 

October 

LT 

February 

May 

August 

November 

00-04 

4.6 

3.7 

5.6 

3.9 

04-08 

4.2 

5.0 

7.0 

4.5 

08-12 

4.5 

6.0 

7.7 

4.9 

12-16 

4.8 

6.3 

9.0 

6.0 

16-20 

4.7 

4.5 

5.9 

4.4 

20-24 

5.3 

3.9 

4.8 

4.3 

Average 

4.68 

4.90 

6.67 

4.67 

Table  3.  Standard  Deviation  (dB)  of  F  Values 

(Op.  )  at  1  MHz  from  CCIR  Report  322. 
am 


Table  3  gives  values  of  °Fam  at  1  MHz  for  each  of  the  seasons  and  time 
blocks.  It  also  gives  the  average  over  time  of  dpam.  Note  that  there  is 
considerably  more  error  in  CCIR  Report  322  in  the  June,  July,  August  season. 
Hence  a  slight  increase  in  rms  residual  in  a  model  for  this  season  might  not 
be  as  important  as  the  error  in  the  data  used  to  generate  the  model. 

Accuracy  of  Zacharisen  and  Jonas  Nodal 

In  the  previous  section.  Table  2  gave  the  rms  error  for  the  Zacharisen 

and  Jones  model.  Table  4  given  below  gives  the  average  rms  error  for  the 

Zacharisen  and  Jones  model  along  with  the  average  error  for  the  CCIR  322 

data.  Examination  of  the  table  shows  the  error  in  the  Zacharisen  and  Jones 

model  is  well  below  that  of  the  CCIR  32  2  data.  Because  the  bias  in  the 

Zacharisen  and  Jones  model  is  nearly  zero,  its  rms  error  is  approximately 

equal  to  the  standard  deviation  of  the  residuals.  Assuming  that  the  rms 

residual  of  F„  of  this  model  and  the  standard  deviation  of  F,_  values  used  to 
am  din 

generate  the  model  are  independent,  then  the  total  error  in  the  model  is  given 
by  the  square  root  of  the  sums  of  the  squares  of  the  two  errors.  This  is  also 
given  in  Table  4. 
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Season 

December 

January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

Average 

rms  residual  (dB) 

1  .62 

1.45 

2.02 

1 .49 

Average 

Op  at  1  MHz 

am 

4.68 

4.90 

6.67 

4.67 

Combined 

4.95 

5.11 

6.97 

4.90 

Table  4. 

Comparison  of  the 

Accuracy 

of  the  Zacharisen  and  Jones 

Atmospheric  Noise  Model  and  the  standard  deviation  of  F 

am 

values  used  to  generate  the  model. 


Accuracy  of  Mlnicoaputar  Models 

Table  5  compares  the  accuracy  of  the  several  minicomputer  models  and  that 

of  the  standard  deviation  o„  used  to  generate  the  models.  The  table  gives 

ram 

the  number  of  coefficients  per  time  block  and  the  longitude  and  latitude 

maximum  harmonic  numbers  J  and  K,  respectively.  As  can  be  seen  from  the 

table,  a  large  number  of  coefficients  (i.e.,  192  of  110)  results  in  a  rms 

error  much  lower  than  that  of  the  CCIR  322  data.  Whereas,  a  small  number  of 

coefficients  (i.e.,  50  and  .35)  results  in  a  rms  error  approaching  that  of  the 

CCIR  322  F  values.  Between  these  two  extremes  are  the  models  with  77  and  55 
am 

coeff icients,  respectively.  The  model  with  55  coefficients  appears  to  be  a 
good  model  when  a  small  number  of  coefficients  is  desired.  Note  that  the 
model  for  192  coefficients  per  time  bloch  is  more  accurate  than  the  Zacharisen 
and  Jones  model,  and  that  the  model  using  110  coefficients  is  more  accurate 
than  the  Zacharisen  and  Jones  model  except  during  the  September,  October, 
November  season. 
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Table  5.  Accuracy  of  Minicomputer  Atmospheric 
Noise  Models. 


INTERPOLATION  KITH  TIKE 


Because  there  is  a  separate  numerical  map  of  F  for  each  of  the  ITT  hours 

*  am 

h  »  0000,  0400,  0800,  1200,  1600,  2000,  it  is  often  necessary  to  interpolate 
between  these  hours  to  get  Fao  at  a  particular  receive  site. 


An  alternate  approach  to  the  standard  linear  Interpolation,  would  be  to 
use  a  finite  Fourier  series  for  Interpolation  In  much  the  same  way  that  the  LT 
data  were  converted  to  UT  data  points.  The  main  difference  here  is  that  the 
formula  for  the  hour  angle  given  in  equation  (11)  is  replaced  by 


tt  -  15  i  ,  1  =  0,  4,  6,  20.  (25) 


By  Fourier  analysis  of  F*h'  for  these  six  hours,  a  continuous  diurnal 

am 

representation  is  given  by 
2 

Y(t  )  -  a  +  I  (a. cos  jt  +b  sin  jt  }  +  a  cos  3  t  (26) 

1  °  j  -  1  3  1  3  13  1 

where  the  Fourier  coefficients  are  given  by 


1  r  w  (h) 
a  -  t  2  F 
o  6  am 
n 


1  „  (h) 

a.  “  t-  1  F  cos  j  t. 

j  3  am  i 


1  „  (h) 

b  «  —  £  F  sin  j  t 

j  3  .  am  J  i 


a,-;*:  (-i,h-1  F.  (h) 


“3  6 


h  -  1 


j-1,2 


(27) 


Here  l  notes  a  sum  over  h  «  0000,  0400,  0800,  1200,  1600,  and  2000.  Having 
h 

obtained  the  coefficients  In  (27),  then  (26)  can  be  used  for  the  hours  0«  4  tt 
<  360*  or  0  <  i  <  24. 
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MODS.  OSS  STB  OTHER  MOISE  SOURCES 


In  the  HF  band  where  these  models  are  Intended  to  be  used,  there  are 
three  principal  types  of  environmental  noise: 

(1)  atmospheric  from  thunderstorms  and  other  types  of  disturbed  weather 

(2)  cosmic  (galactic  noise) 

( 3 )  man-made 

To  assess  the  effect  of  radio  noise  on  a  receiver  system  it  is  necessary  to 
calculate  the  combined  noise  power  available  from  an  equivalent  loss  free 
antenna. 

Atmospheric 

Atmospheric  noise  as  calculated  by  these  models  is  usually  the  dominant 
type  of  noise  at  HF,  especially  in  the  lower  portion  of  the  band.  It  is  also 
the  most  erratic  of  the  three  major  types  of  noise.  It  may  change  over  wide 
limits  as  a  function  of  location,  frequency,  bandwidth,  time  of  day,  season, 
solar  activity,  and  azimuthal  direction.  Although  in  the  presence  of  local 
storms  atmospheric  noise  may  be  an  important  factor  at  almost  any  frequency: 
it  is  the  ability  of  noise  from  distant  thunderstorms  to  propagate  over  long 
distances  that  makes  it  so  important  at  HF  and  below.  Because  ionospheric 
absorption  is  high  during  the  daytime,  the  contribution  from  distant  sources 
is  reduced,  and  local  sources  become  important.  Because  of  the  strength  of 
propagated  noise  from  distant  storms  at  night,  the  diurnal  maximum  occurs  at 
night  -  even  for  locations  in  the  earth's  major  source  regions. 

For  the  purposes  of  combining  the  component  noise  into  one,  denote 
atmospheric  noise  at  a  particular  frequency  by  Na  and  its  variability  by  Dua 
and  D^a- 

Cosmic 

Cosmic  noise  is  thermal  noise  that  originates  outside  the  earth  and  its 
atmosphere.  It  includes  both  solar  radio  noise  and  noise  from  interstellar 

space. 
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The  low  frequency  cutoff  of  cosmic  noise  Is  determined  by  the  peak 
density  of  the  Ionosphere;  radio  waves  of  frequency  less  than  the  critical 
frequency  of  the  F-region  cannot  penetrate  to  the  earth.  Cosmic  noise  starts 
to  become  an  important  factor  near  the  frequencies  where  propagated  atmos¬ 
pheric  noise  drops  due  to  lack  of  ionospheric  support  and  absorption.  Because 
the  existence  of  cosmic  noise  at  the  ground  depends  upon  the  critical  fre¬ 
quency  of  the  F-reglon,  it  is  necessary  to  estimate  the  critical  frequency  in 
order  to  determine  whether  cosmic  noise  should  be  considered.  On  the  other 
hand,  it  is  not  desirable  to  Insert  a  major  portion  of  an  HF  skywave  predic¬ 
tion  program  in  a  minicomputer  to  see  whether  galactic  noise  should  be  con¬ 
sidered.  Because  cosmic  noise  is  rarely  the  dominant  type  of  noise  for  an  HF 
system  in  the  ambient  environment  except  at  the  highest  frequencies  and  at 
very  quiet  receiving  sites,  it  suffices  to  determine  the  maximum  expected 
median  fQF2  for  low,  moderate,  and  high  sunspot  levels  for  each  season.  This 
was  accomplished  by  examining  contour  plots  of  fQF2  for  sunspot  numbers,  10, 
110,  and  160. 14-17  The  maximum  contour  value  for  each  month  was  determined 
for  each  sunspot  level.  Then  the  average  of  the  three  monthly  values  was 
chosen  as  the  value  to  be  used  for  the  maximum  expected  median  fQF2 .  The 
values  so  determined  are  presented  in  table  6.  These  values  are  to  be 
interpolated  for  the  actual  sunspot  number,  and  cosmic  noise  is  then  to  be 
considered  if  the  operating  frequency  is  above  this  interpolated  frequency. 


Sunspot 

Number 

December 

January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

10 

11.3 

12.3 

10.  P 

12.0 

110 

15.0 

15.3 

13.7 

16.7 

160 

16.0 

16.0 

14.0 

17.3 

Table  6. 

Maximum  Median 

MUF ( ZERO) F2 

At  times  when  it  is  possible  to  see  the  galaxy  at  the  operating  fre¬ 
quency,  the  median  galactic  noise  is  given  for  a  short  vertical  electric 
dipole  above  a  perfectly  conducting  ground  screen  by 


*  52.0  -  23.0  l°g10  f 


(28) 
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=  expected  median  value  of  the  galactic  noise  in  a  1  Hz  band 
— dB  >  1  watt 

f  =  operating  frequency  —  MHz .  ^ 

The  variability  of  galactic  noise  !DUg  and  Dig^  can  **  taken  to  2  dB 
about  the  median. 

Han-made  Noise 

At  certain  receiving  locations,  unintended  man-made  radio  noise  may  be 
the  predominant  external  noise  with  which  the  communication  system  must 
compete.  It  may  arise  from  any  number  of  sources,  such  as  power  lines, 
industrial  machinery.  Ignition  systems,  etc.,  and  thus  it  may  have  wide 
geographic  and  short  term  variations.  Because  it  propagates  primarily  by  the 
ground  wave  and  along  power  lines,  its  range  of  effects  is  limited  and  its 
average  value  depends  primarily  on  the  location  of  the  receiving  site  relative 
to  the  sources. 

A  major  advance  in  environmental  noise  level  documentation  was  made  when 
workers  at  the  U.  S.  Department  of  Commerce  Institute  for  Telecommunications 
Sciences  (ITS)  used  the  same  measurement  system  to  obtain  noise  data  in  the 
band  250  kHz  to  250  MHz  with  a  short  vertical  antenna  near  ground  at  various 
sites  in  the  United  States. Three  environmental  categories  were  desig¬ 
nated:  rural,  residential,  and  business.  Rural  areas  were  defined  as  loca¬ 

tions  where  land  usage  is  primarily  for  agricultural  or  similar  pursuits,  and 
dwellings  are  no  more  than  one  every  five  acres.  Residential  areas  (urban  or 
suburban)  were  defined  as  any  area  used  predominantly  for  single  or  multiple 
fully  dwellings  with  a  density  of  at  least  two  family  units  per  acre  and  no 
large  or  busy  highways.  A  business  area  was  defined  as  any  area  where  the 
predominant  usage  throughout  the  area  is  for  any  type  of  bsiness  (shopping 
centers,  main  streets  or  highways  lined  with  various  business  enterprises, 
etc.).  Further  limited  measurements  relate  to  parks  and  university  campuses 
(which  generally  have  less  vehicular  traffic,  electrical  equipment  and  other 
types  of  noise  sources  than  business  or  residential  areas)  and  relate  to 
interstate  highways  (where  the  main  source  is  ignition  noise) . 


The  results  of  the  measurements  were  analyzed  statistically  and  least- 
squares  fits  for  F  were  obtained.  In  all  cases  results  were  consistent  with 
a  linear  variation  of  the  median  value,  tMm,  with  frequency  f  of  the  form 

Mm  -  c  -  d  log10  f  (29) 

with  f  expressed  in  MHz,  and  c  and  d  taking  the  values  given  in  table  7. 

18 

Using  the  data  given  by  Spaulding  and  Disney,  Hagn  and  Sailors  found  that 
the  time  and  location  variability  of  man-made  noise  could  be  represented  by  a 
simple  Gaussian  model.  The  standard  deviations  for  this  model  are  given  in 
table  8. 


Environmental  Category 

c 

d 

Business 

76.8 

27.7 

Interstate  Highways 

73.0 

27.7 

Residential 

72.5 

27.7 

Parks  and  University  Campuses 

69.3 

27.7 

Rural 

67.2 

27.7 

Table  7.  Values  of  the  constant  c  and  d  in  the  Equation 

for  F  for  Man-Made  Noise 

am 


Freq. 

Business 

Residential 

Rural 

(MHz) 

Pam 

°NM 

<dB) 

F 

am 

°NM 

(dB> 

F 

am 

°NM 

(dB) 

1.00 

76.8 

6.3 

72.5 

6.5 

67.2 

9.5 

2.50 

65.8 

12.4 

61.5 

10.4 

56.2 

10.2 

5.00 

57.4 

9.2 

53.1 

8.4 

47.8 

9.4 

10.00 

49.1 

7.7 

44.8 

6.1 

39.5 

6.8 

20.00 

40.8 

8.7 

36.5 

8.3 

31.2 

6.9 

48.00 

30.2 

11.1 

25.9 

8.8 

20.6 

4.5 

Table  8. 

Summary 

of  Man-Made 

Noise  Model 

Parameters. 
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combination  of  Noise 


Vfhen  the  noise  received  is  due  to  several  different  sources,  the  com¬ 
ponent  noise  will  combine  into  one,  talcing  the  different  phases  into  account. 
In  this  case,  these  phases  are  independent  and  continuously  variable.  Hence, 
the  noise  power  in  watts  can  be  summed.  The  resultant  median  noise  is  given 
by 

N  /10  N  /10  N  /10 

N  -  10  iog1()  (10  a  +  10  9  +  10  “  )  (30) 

where  Na  is  the  atmospheric  noise  in  a  1  Hz  band  --  dB  above  1  watt.  Figure 
20  shows  the  quadratic  addition  chart  for  two  signals.  If  two  noise  sources 
are  within  3  dB  of  each  other,  then  an  additional  1.76  dB  of  noise  is 
received.  On  the  other  hand,  if  the  two  sources  are  within  only  6.02  dB  of 
each  other,  the  additional  noise  is  0.976  dB.  Finally,  if  the  two  sources  are 
separated  by  12.04  dB,  the  additional  noise  due  to  a  second  source  would  be 
0.26  dB. 

The  upper  and  lower  decile  differences  from  the  median  noise  are  found 
also  by  summing  the  upper  and  lower  deciles  of  the  noise  in  watts,  converting 
the  deciles  back  to  decibels,  and  subtracting  the  median  noise.  The  resulting 
expressions  are: 


D  » 


HU 


(N  *  D  )/10 

1 10  log1Q(10  a  u*  +  10 


(N  +  D  )/10 
g  ug 


(N  +  1.280  )/10) 

10  m  NM 


(31) 


Hi 


|101og1(J(10 


(V 


D,  )/10  (N  +  D,  ) 1 0 
la  +  10  9  19  ♦ 
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where  Dua  an<j  D  are  the  Upper  an<j  lower  deciles  of  atmospheric  noise.  The 
results  of  equations  (31)  and  (32)  can  be  used  to  determine  circuit 
reliability. 
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dB  DIFFERENCE  BETWEEN  TWO  SIGNALS 


Figure  20.  Quadratic  addition  chart  for  two  signals  in  dB 
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CONCLUSIONS 


Practical  but  simplified  atmospheric  noise  models  were  developed  for  use 
in  estimating  system  signal-to-noise  ratios  in  micro-mini-computer  based  HF 
propagation  prediction  systems-  The  numerical  mapping  techniques  of 
Zacharisen  and  Jones  were  used  to  represent  worldwide  atmospheric  noise  at  1 
MHz.  Each  four-hour  time  block  was  represented  separately  rather  than  for  a 
twenty-four  hour  day.  This  resulted  in  24  numerical  maps:  one  for  each  of 
six  4-hour  time  blocks  (in  universal  time,  UT)  and  each  of  four  3-month 
periods.  This  method  provided  a  numerical  representation  of  atmospheric  noise 
with  a  minimum  reduction  (by  one-fifth)  from  the  number  of  coefficients  used 
by  Zacharisen  and  Jones  in  computer  storage  at  one  time.  The  remaining 
coefficients  would  be  stored  in  the  computer's  cassette  tape  unit.  The 
simpler  expression  for  the  median  noise  levels  (F  )  results  in  significant 
savings  in  computer  code  and  time.  Several  versions  were  produced-each  with  a 
different  number  of  coefficients  (number  of  harmonics)  in  the  Fourier 
representation  in  each  time  block.  The  models  produced  make  use  of  the 
frequency  dependence  model  of  earlier  representations  of  CCIR  Report  322. 

The  accuracy  of  these  model  versions  was  determined  as  a  function  of  the 

number  of  harmonics  in  the  representation.  The  average  (bias)  and  rms 

residual  was  determined  between  values  of  F _  used  to  develop  the  model  and 

am 

the  corresponding  values  computed  from  each  numerical  model.  In  all  cases  the 
bias  was  nearly  zero.  The  significance  of  the  rms  residual  was  determined  by 
comparing  the  value  of  each  version  to  the  published  accuracy  of  the  CCIR  322 
noise  data  used  to  generate  the  models. 

The  average  rms  error  for  the  Zacharisen  and  Jones  model  for  all  hours 
and  seasons  was  1.65  dB.  Its  worst  season:  June,  July,  and  August:  had  an  rms 
error  of  2.02  dB.  The  corresponding  rms  errors  for  the  CCIR  data  were  5.48  dB 
and  6.97  dB,  respectively.  The  Zacharisen  and  Jones  model  rms  error  is  much 
smaller  than  the  CCIR  322  data. 

The  average  rms  error  for  all  times  and  seasons  for  the  model  using  192 
coefficients  per  time  block  was  1.25  dB.  Its  worst  season  was  the  June,  July, 
August  season,  which  had  an  rms  error  of  1.4  dB.  These  errors  were  well  below 
that  of  the  CCIR  322  noise  data  and  are  even  less  than  that  of  the  Zacharisen 
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and  Jones  model.  The  most  suitable  application  for  this  version  would  be  in 
Advanced  Prophet  applications  or  in  a  large  scale  HF  prediction  model. 


The  average  rms  error  for  all  times  and  seasons  for  the  versions  with  110 
'-oef f iclents  was  1.60  dB.  For  the  June,  July,  August  season  the  error  was 
1.76  dB.  These  error  figures  were  well  below  that  of  the  CCIR  322  data  and 
were  slightly  less  than  that  of  the  Zacharisen  and  Jones  model.  The  reason 
this  model  had  better  accuracy  than  the  Zacharisen  and  Jones  model  was  that  it 
retained  more  time  information.  The  most  suitable  application  for  this 
version  would  be  use  in  Advanced  Prophet  applications. 

The  two  versions  with  77  and  55  coefficients  each  had  an  average  rms 
error  for  all  times  and  seasons  of  2.18  dB  and  2.94  dB,  respectively.  During 
June,  July,  August,  these  two  versions  had  a  rms  error  of  2.8  dB  and  3.57  dB, 
respectively.  These  error  figures  were  still  reasonably  below  that  of  the 
CCIR  322  noise  data.  These  models  would  be  suitable  for  both  Advanced  Prophet 
applications  and  for  applications  using  Tektronix  4052  type  mini-computers. 

Finally,  the  two  versions  with  50  and  35  coefficients  each  had  an  average 
rms  error  for  all  times  and  seasons  of  3.64  dB  and  4.53  dB,  respectively. 
During  June,  July,  August,  these  two  versions  had  a  rms  error  of  5.25  dB  and 
5.63  dB,  respectively.  These  error  figures  have  begun  to  approach  the  average 
rms  error  of  5.23  dB  for  the  CCIR  322  noise  data.  These  versions  would  be 
most  suitable  for  applications  for  which  the  memory  requirements  for  the  mini¬ 
computer  demand  a  sacrifice  in  the  accuracy  of  the  model. 
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RECOMMENDATIONS 


As  a  result  of  this  study,  the  following  recommendations  are  made: 

1.  Fourier  series  interpolation  should  be  used  to  obtain  F^  at  1  MHz 
for  an  actual  hour  of  operation; 

2.  When  Fourier  series  interpolation  is  not  practical,  then  linear 
interpolation  should  be  used; 

3.  Because  other  noise  sources  are  present  at  KF,  these  noise  values 
should  also  be  determined; 

4.  Galactic  noise  should  only  be  calculated  for  frequencies  above  the 
maximum  expected  fQF2  given  herein  as  a  function  of  sunspot  number  and  season, - 

5.  The  formula  provided  should  be  used  for  combining  the  values  of  man¬ 
made,  galactic,  and  atmospheric  noise; 

6.  The  reliability  of  the  combined  1  MHz  models  and  the  frequency 
dependence  models  should  be  determined; 

7.  This  reliability  should  be  determined  as  a  function  of  frequency, 
time  of  day,  geographic  location  as  well  as  the  maximum  longitude  and  latitude 
harmonic  numbers  of  the  1  MHz  numerical  maps; 

8.  The  accuracy  of  the  complete  Lucas  and  Harper  Model2  should  be 
determined  with  the  same  data; 

9.  This  reliability  analysis  should  make  use  of  data  measured  after  the 
data  contributing  to  CCIR  Report  322  were  measured  as  well  as  for  data  used  to 
develop  CCIR  Report  322. 
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APPENDIX  A 


WORLD  CONTOUR  MAPS  OF  1  MHZ  ATMOSPHERIC  NOISE 

UT  contour  maps  of  Fam  In  dB  above  KTQb  are  particularly  useful  In  making 
choices  between  the  various  1  MHz  atmospheric  noise  models*  The  absence  of 
contours  of  high  level  or  the  lack  of  geographical  detail  is  particularly 
significant. 

The  June,  July,  August  season  and  December,  January,  February  season 
were  chosen  to  illustrate  the  effects  of  reducing  the  number  of  coefficients. 
Contour  maps  were  drawn  for  the  UT  hours  0000  and  1200.  A  contour  map  was 
drawn  for  each  harmonic  number  combination  of  J=10  and  K=9;  of  J=*5  and  K=9,  6, 
4;  and  J=2  and  K-9 ,  7.  A  contour  plot  was  also  drawn  for  each  hour  and  season 
for  the  Zacharlsen  and  Jones  model.  Figures  13-19  in  the  main  body  of  the 
report  are  for  the  hour  =  0000  during  June,  July,  and  August.  The  remaining 
contour  maps  are  presented  here  in  figures  A1-A21. 

Figures  A1  through  A7  are  for  December,  January,  February  at  0000  hours 
UT.  Figure  A1  is  for  the  Zacharisen  and  Jones  model.  The  remaining  six 
figures  are  for  the  six  versions  of  the  1  MHz  atmospheric  noise  model 
presented  in  the  main  text. 

Figures  A8  through  A14  are  for  December,  January,  February  at  1200  hours 
JT.  Figure  A8  is  for  the  Zacharisen  and  Jones  model,  and  Figures  A9-A14  are 
for  the  six  model  versions. 

Figures  A15  through  A21  are  for  June,  July,  August  at  1200  hours  OT. 
Figure  A15  is  for  the  Zacharisen  and  Jones  model,  and  Figures  A16-A21  are  for 
the  six  model  versions  being  compared. 

The  same  type  of  variation  from  contour  map  to  contour  map  reported  on 
Figures  13-14  also  occurs  here.  For  the  map  versions  with  J=10,  K*9,  there  is 

slightly  more  detail  in  the  equatorial  and  polar  regions  than  in  the 

Zacharisen  and  Jones  model.  The  model  for  J=5,  K=9  is  very  much  like  that  for 

the  Zacharisen  and  Jones  model  except  there  is  a  little  more  detail  in  the 
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polar  regions.  As  the  number  of  harmonics  Is  reduced  further,  the  high  level 
contours  begin  to  disappear.  In  the  June,  July,  and  August  season,  this 
change  occurs  even  with  higher  order  harmonics.  For  the  J  =2,  K  =  9  and  for 
the  J  =  2,  K  =  6  cases,  there  is  a  loss  of  detail  over  the  whole  contour  plot 
from  the  versions  with  larger  harmonics. 
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APPENDIX  B 


REPRESENTATION  OP  FREQUENCY  DEPENDENCE  AND  VARIABILITY 
OP  ATMOSPHERIC  1*31  SE 


PURPOSE 

To  determine  the  value  of  for  atmospheric  noise  for  a  given  time  and 
frequency  and  its  variability* 

NUMERICAL  METHOD 

The  frequency  Dependence  curves  given  in  CCIR  Report  322  were  fitted  by 
Lucas  and  Harper  using  least  squares.2  The  result  is  a  power  series  of  the 
form 


Y(X,N)  -  A^ (N)  +  A2(N)  X  +A3(N>  X2  +  ... 
+  A?(N)  X6 


(B-1) 


where 


VH)  *  bi,1  **  1  “  1* 


(B-2) 


N  “  amplitude  of  the  hourly  median  1  MHz  atmospheric  noise  (dB>kTQb) ,  and 


X 


logio  £ 

8x2  -11 


(B-3) 


76 


where  f  is  the  operating  frequency  in  MHz.  Equation  (B-1 )  is  subject  to  the 
constraint  ¥(-3/4,11)  •  N  (i.e.,  the  frequency  variation  of  the  1  MHz  value 
must  be  zero) . 


The  CC1R  Report  322  curves  representing  the  decile  values  of  the 
distributions  were  generated  by  a  least  square  fit  by  Lucas  and  Harper  to  the 
function 

2  3  4 

Y(X)  «  A1  +  A2  X  +  A3  X  +  A4  X  +  Aj  X  (B-4) 

where  X  *  10  l°g10  f  end  f  is  the  operating  frequency  in  MHz. 


STRUCTURE 

SUBROUTINE  GENFAM(RLAT,  NTB,  FREQ,  AHOS,  ATNO,  DU,  DL,  GU,  SIGU,  SIGL,  SIGM) 


Formal  Parameters 

RLAT 

NTB 

FREQ 

ANOS 


Geographic  latitude  in  radians  (+N,-S) 

Time  band  index 
Frequency  in  MHz 

1  MHz  atmospheric  noise  for  time  of  interest 


ATNO,  DU,  DL 


Atmospheric  noiBe  and  deciles  at  desifed 
freque  cy 


SIGU,  SIGL,  SIGM 


Standard  deviation  of  deciles  and  mean  noise 


Related  Algorithms 

An  algorithm  called  XINTER  provides  the  interpolation  time  band  Indices 
for  use  with  GEN F AM  and  the  linear  interpolation  fraction  for  a  given  UT. 
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SUBROUTINE  XINTER  (TIMER,  I  TIME,  NTB,  ABCT ) 


Formal  Parameters 

TIMER  Time  —  either  local  or  universal 

I TIME  Time  type  index  —  1  for  LT,  2  for  UT 

NTB  First  time  band  index  for  interpolation 

NTX  Second  time  band  index  for  interpolation 

ABCT  Fraction  for  interpolation 

Usage 


For  use  with  GENFAM  it  is  necessary  to  determine  the  UT  corresponding  to 
the  UT  of  interest  at  the  receive  site.  Then  XINTER  is  called.  Next  two 
calls  are  made  to  GENFAM  —  one  with  NTB,  one  with  NTX  (1  i.e.,  with  NTB 
atmospheric  noise  is  ATNO;  with  NTX  atmospt’x  :  noise  is  ATNX).  Then  the 
atmospheric  noise  at  a  hour  LT  corresponding  OT  is 


ATMNO  -  ATNO  +  (ATNX  -  ATNO)  *  ABCT  ( B- 5 ) 

In  the  same  manner  the  upper  and  l<x»er  deciles  are  given  by 

DUA  -  DU  +  (DUX  -  DU)  *  ABCT 

(B-6) 

and 

DLA  *  Dli  t  (DLX  -  DL)  «  ABCT  (8-7) 
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SUEPOUTINE  GENFAM  (R|.AT,NTP*FRE0»AMnS,ATH0,PU,t'|.  ,SlPU,SIGL»S!Gt'  ) 

HUA03100 

c 

c 

CALCULATE®  THF  FREQUENCY  DEPENDENCE  OF  the  MFrtiN 

1  mu7  ATMOSPHERIC  noise  and  its  associated  decile  values  and 

c 

c 

STANDARD  DEVIATION  OF  THE  MEDIAN  AND  UPPER  aK'T  LOWE#  DfCTCFS 

CF  THE  ATMOSPHERIC  NOISE. 

c 

c 

RLAT*PECGRAPHIC  LATITUDE  IN  RAD  I  AMS ( +N.-S » 

NTS«TIME  PAND  INDEX 

c 

c 

ereo*fpeouency  IN  MH7 
anos»i  yhy  atmospheric  noise 

c 

c 

ATNO,rL’,PL*ATMOSPHERIC  NOISE  AND  DECILES  AT  DESIRED  FREQ 
MGU.SIGL.SIGM-STD.  DEVIATIONS  OF  DECILES  AND  MEAN  NOISE 

IHRLltrT  REAL*6  U-H.0-7) 

DIMENSION  V ( 5  ) 

HUA03101 

COMMON  /ONE /  Dl'D(5, 12, 5)»EAH<  14,12) 

IF  (FPFO  .GT.  30)  GO  TC  6 

Nil  AOS  105 

F20-»LOG10(20.0) 

F«AL0G10'FREC) 

HUA03106 

HUA03107 

itr*ntr 

IF  (PLAT. LT. o.O >  ITP-ITR+6 

MUA03108 

HUA03109 

1 

X  — 0.75 

P7«0.0 

HUA03110 

HUA031I1 

PV. o.O 

CO  2  I  *  t  »  7 

MUA03112 

HUA03113 

2 

P  7*  X  *  P7  +  F  A  H ( I , ITR ) 

PY«X«Py+FAM ( 1+7, 1TB) 

NUA03I14 

HU A 031 15 

IF  (X.NE.-0.75)  GO  TC  3  ” 

C7»aN0S*(2.0-PZ)-PX 

HUA03116 

HU  A 031 1  7 

X«(P.0*2.0**E-H.o>/4.0 

GO  TO  l 

HUA0311A 

HU  A 031 1 9 

3 

ATNQ*C?*P?+PX 

IF(F.GT.F20)  F-F20 

HUA03120 

HUA03121 

DO  5  1*1,5 

Y-0.0 

HUA0312?  ~ 
HUA03123 

IF( I.E0.5.AND.F.GT.1.0I  F«1.0 

DO  4  J*  1 , 5 

^U A 03 124 
HUA03125 

4 

5 

Y-F,Y*0I.C  <  J,  ITP.  t  ) 

V(  J)*Y 

^UAo3i26 

HUA03127 

CU*V  (  1  ) 

TL*v(2l 

SIGU»V(3)  “ 

S I GL* V ( 4 ) 

HUA03130 

HUA0313I 

SlCM.V(f)  ~  “ 

RFTURN 

“HLT0TI35 - 

HLA03133 

6 

atKo  -  -10  "  '  - - - - 

DL  *  0 

“^uxrrrnra — 

HUA03135 

DL  *  0 

SIGU-  0 

HLA0313F 
HUA03137  ‘ 

§,GL  .  0  -  ■  —  -  -  —  - 
SIGN  «  0 

“VT3T03I3P 

HUA03139 

—  rfYupn  “  ‘  " 

EN[ 

“  ^L)  A03 1 40  ' 

HUA03141 

atmospheric  noiaa 
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_ SUBROUTINE  V  INTER < TIMER, IT IMF, NTP#NTX.APCT) _ 

c  provides  interpolation  indices  ant  fraction  for  sane 

C  TINER«TIHE--EITHER  local  or  universal  tine 

C  ITIME-TINE  TPYE  INDEX  — l  FOR  LT.  2  FOR  UT 

c _ NTB»FIRST  tine  rand  index _ 

C  NT*«SECCNT  TIME  RA^D  INDEX 

C  APCT»FRACTION  FQR  INTERPOLATION _ 

IMPLICIT  RE AL*6  (A-H.O-7) 

_ CC»TIMFR _ 

NTReCC/4. 0*1.0 

_ I F I NT6  .GT.  6)  NTR«6 _ 

TM«4*NTP-2*!TINE 

_ AgCT«APS<CC-TM»/4.0 _ 

NTVaNTB-t 

_ IF  ICC  .GT.  TMI  NTXeNTB*  1 _ 

IFINTX  .GT.  6)  NTX« l 

_ IFINTX  .LT.  1)  NTX"6 _ 

RETURN 

_ END _ 


Figure  B2  Listing  of  XXNTXR,  an  algorithm  for  determining  linear 
Interpolation  Indices  and  fraction  for  same 
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APPENDIX  C 

FORTRAN  SUBROUTINE  OP  THE  1  MHZ 
ATMOSPHERIC  NOISE  MODEL 


PURPOSE 

To  calculate  the  value  of  for  atmospheric  noise  for  a  given  time  at  1 

MHz. 


NUMERICAL  METHOD 

The  two-dimensional  geographic  representation  of  F  for  the  hours  0000, 

am 

0400,  0800,  1200,  1600,  and  2000  for  each  season  is  given  by 

K  (h) 

X  (A,  0)  »  l  p  G  (X) 

h  k»0  o,k  o,k 

J  ,  K  (h) 

+  E  L  (  E  P  G  (X)  cos  j©  (C-1) 

j-1  k-0  j,k  j,k 

K  (h) 

+  (  E  q  G  (X))  sin  j© 

K-0  j,k  j,k 

where  maximum  longitudinal  harmonic  number,  k,  is  9  for  0<  j  <  7  and  is  6  for 
j  *  8,  9,  10,  and  where  J  is  the  longitudinal  maximum  harmonic.  The  latitude 
functions  X)  are  defined  by 

Gj,k  (A)  “  ®ink  A  cos^  A  (C-2) 
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Here  X  denotes  latitude,  -90*  <  X  <90*,  and  6  denotes  longitude  measured 
eastward  from  0*,  0*  <  6  <  360°. 

The  coefficients  p*h*  and  q^h'  were  found  using  the  numerical  mapping 

Jr*  2 

techniques  developed  by  Zacharisen  and  Jones*  The  procedure  is  described  in 
the  main  body  of  the  report.  The  actual  number  of  coefficients  used  depends 
on  the  memory  space  available  and  the  application. 


STRUCTURE 

The  model  consists  of  two  subroutines  -  EVALGK  and  EVAL4  listed  in 
Figures  C1-C2.  The  first  evaluates  the  G^k(X)  functions  and  the  resulting 
summation  of  the  product  of  G^k(X)  and  Pj/k*h*  and  ^  in  equation  (C-1) 
at  a  fixed  latitude  (e.g.,  JL  p*h*  G.  .  (X)).  The  latter  subroutine 
completes  the  evaluation  of  Fu  in  equation  (C-1)  at  a  fixed  longitude  using 
the  Fourier  coefficients  generated  in  EVALGK. 


SUBROUTINE  EVALGK (FLAT,  MAXG,  NH) 


Box 


il  Par* 


iters 


FIAT  Geographic  latitude  in  degrees  (+N,  -S) 

MAXG  Maximum  latitudinal  harmonic  K+1 

MM  Maximum  longitudinal  harmonic  J 
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SUBROUTINE  VOUdUM,  n,  OMEGA) 


Formal  Parameters 

FLON  Geographic  longitude  In  degrees  eastward  from  0° 

NH  Maximvm  longitudinal  harmonic  J 

OMEGA  Estimated  value  of  for  given  latitude,  longitude, 

and  time  of  day 

OOMKH  /A/  D(21 , 10)  ,DK(21 )  ,DR 

formal  Parameters 

D(21 , 10 )  Set  of  numerical  coefficients  representing  Fam 

DK (21)  Resultant  Fourier  coefficients  calculated  in  EVALGK 

DS  Degrees  to  radians  constant  (0.01745329252) 


RELATED  ALGORITHMS 


An  algorithm  called  XINTER  provides  the  interpolated  time  band  indices 
for  use  with  these  two  subroutines.  Using  XINTER  it  is  possible  to  linearly 
interpolate  between  the  VS  hours  0000,  0400,  0800,  1200,  1600,  and  2000. 

I  TIME  should  be  set  equal  to  2.  XINTER  and  its  usage  are  described  in 
appendix  B. 

ADDITIONAL  OOHNEMTS 


The  listings  given  for  EVALGK  and  EVAL4  are  in  the  Harris  FORTRAN  66 
language.  The  Harris  FORTRAN  66  language  is  a  super  set  of  USA  FORTRAN  X3.9 
1966. 
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Thr  listings  given  here  are  for  J**10  and  K«9.  For  other  versions  of  the 
numerical  map  it  may  be  necessary  to  make  minor  changes  to  the  subroutines. 
The  dimensioning  in  the  common  block  A  for  D  should  be  (2*NH+1 ,  K+1 )  and  OK 
should  be  (2*NH+1).  In  EVAL4  the  dimensioning  for  FCOS  and  FS1N  should  be 
<2*NH+1). 
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_ _SUeROUT[*'F  FVALGK(FlAT*MAXG»NH) 

C  EVALUATES  THE  GK  FUNCTION. 

IMPLICIT  RE  AL*6  (A-F.0-7) 

CON HON  /A/  D(2|,I0I,DK(2I I.CR 
NF«2*NH+1 
HAT  =  HAXP 
TL AT  ^__FLA_T  *  OR 
’  TSIN  •$  I  N  f  TL  AT ) 

TCCS  ■  COS ( TL  AT ) 

CO  1  I  =  r.NF 
JTERH«i/2 
SUV  ■  o.O 

_ I  F(I  .GE»  1 6  )  HAXG«=7 

CO  2  K*l,HAXG 

KK.K-1 

SUN  ■SUN  ♦  C(I.K)  *  (TSIN  **KK)  *  (TcCS  **  JTERH ) 

2 _ CONTINUE 

1  FK(!)  -  SLiy 

_ HAX£«HAV 

RFTURN 

ENT 


Figure  Cl  Listing  of  EVALGK,  an  algorithm  for  obtaining  1  MHz 
atmospheric  noise 
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SUBROUTINE  EVAL4  {Fl.CN,NH,OFFPA> 

EVALUATES  FOURIER  COEFF  10 IFN  T$  FCR  A  Fixrn  L*T  I  TUtE 
IMPLICIT  RE  *L*6  (A-H.0-7) 

COMMON  /A/  r«2U10),tK(2!  >.ER 
rif'FNSICN'  FCCS(2l  »,FS  IK  (21  > 

SUN  *  D KM) 

RAT  *  HR  *  FLOK1 
CF  *  COSfRAC) 

Fcrs< 1 )  *  CF 
SF  ■  SIN(RAD) 

FS  I  N  < i I  *  SF 
rc  i  i*2, k>- 

Frrs»n  *  cf  *  Fcns(  i  -  i i  -  sf  *  fsinu  -  1) 

FSINU)  *  CF  *  FSINU  -  J  )  ♦  SF  *  FCOSII  -  1) 

TO  2  J*1»NU 

SU*  *  SUN  ♦  rK(2*  J)  *  FSINIJ)  ♦  rK(2  *  .T  +  1  )  *  FCCS(J) 

ONEGA  *  SIN 

RF7UPN 

FNT 


Figure  C2  Listing  of  EVAL4,  an  algorithm  for  obtaining  1  MHz 
atmospheric  noise 
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APPENDIX  D 


A  COMPUTER  PROGRAM  TOR  READING  AND  PRINTING 
THE  ATMOSPHERIC  NOISE  NUMERICAL  COEFFICIENTS 


Purpose 

To  read  into  a  computer  the  atmospheric  noise  coefficients,  generated  by 
UT  NEQSOL,  and  to  print  these  coefficients  in  tabular  listings. 

Method 


The  number  of  cards  produced  by  UT  NEQSOL  depends  on  the  value  of  K,  the 
maximum  latitudinal  harmonic.  Each  card  contains  an  index  identifying  the 
season  and  hour,  an  index  identifying  the  longitudinal  harmonic,  an  index 
identifying  the  beginning  latitudinal  harmonic  on  the  card,  and  four  Fourier 
coefficients.  It  produces  three  cards  for  the  first  ten  latitudinal  coeffici¬ 
ents  for  J  .LT.  16  and  two  for  J  .GE.  16.  This  results  in  342  cards  per 
season  when  K.-9  and  252  when  K=4  or  6. 

Input 

The  input  consists  of: 

1)  An  alphanumerical  title  card  (columns  1-78) 

2)  A  card  with  J,  K+1,  season  number,  and  a  control  index  to  control  the 
end  of  the  program. 

3)  The  set  of  coefficient  cards. 

Each  subsequent  season  is  placed  after  the  one  preceding  it.  The  variable 
names  and  format  of  the  input  are  given  as  follows: 


1) 

NTITLE 

13A6 

Title  to  be  printed 

2) 

NH 

110 

J 

MAXG 

110 

K+1  (e.g. ,  10  for  K*9) 

NSEAS 

110 

Season  number  (1  through  4) 

NEND 

110 

.GE .10  will  halt  program 
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O DROI 


The  output  looks  like  that  shown  in  Figure  D- 1 .  Across  the  top  is  the 
title.  Then  hours  0000,0400,0800,1200,1600,  and  2000  appear  as  column  head¬ 
ings.  Bach  column  is  for  a  particular  hour  UT.  On  the  left  are  values  of  j 
and  k.  Each  row  is  for  a  combination  of  the  values  of  j  and  k.  The  first  row 
is  for  j~1  and  k=*1 .  Subsequent  rows  are  for  increasing  k  values  up  to  10. 
Then  j  is  incremented  and  k  starts  again  at  1 .  This  process  continues  until 
both  j  and  k  reach  their  maximum  value  (2*J+1  and  K+1). 


Since  UT  NBQSOL  considers  J  to  be  10,  the  printout  provided  by  this  pro¬ 
gram  continues  until  j  equals  10.  Some  of  the  models  described  in  the  main 
body  of  the  report  only  have  a  J  equal  5,  As  all  the  longitudinal  coeffici¬ 
ents  are  orthogonal,  the  longitudinal  coefficients  for  j  values  greater  than 
2*J+1  can  be  ignored  for  that  model. 
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C  READS  Tv-F  NUFER  I  A(_  COEFFICIENTS  FROM  IJT  NECSOL 

~£  AN  t  PRINTS  TNEN  FOR  EACH  SEASON 

c  favimum  value  for  nh  is  io 

C  F  A  *  I  NUN  VALUE  FOR  Na'xg  IS  10 

IFPLICIT  REAL*6<A-H.C-7> 

RFAL*6  NTITLE 

_ r  I  FENS  I  ON  NTITLEI  13)  JFM2U1.0»6) _ 

1  READ  100.  CNTITLEm.IM.l3) 

2  READ  2OO0  NH.MAXG.NSEAS.NENr 
NF»2*NH+1 

NCAROS-57 

lP(HA)fO  ."IE.  M  NCARCS-42 

_ ro  3  i«i,k  _ 

CO  3  J»1 »N CARDS 

READ  4no,NFILE.NHi,NFF,FN(NM.NFF»n.FN(NHl.N'EFM.  1  >, 

1  FN  (NHr.NFF*2,I),FN<NHl,NFF>3,  I  I 

3  CONTINUE 

9  PRINT  500 

_ I  PR  1 NT«U  _ 

CO  700  I-I.NF 
FF»“AXG 

ITTT75F7nrT  fP»? 

CO  300  JM.FF 

IPRINt-lPRINt+1 

IFfIPRIN'T  .GT.  4M  PRINT  500 

- 1  f wrz gt "aw  t  •} pb  r«T- \ nn nrr=i e - 

I F  (  I  PR  I  NT  .EC.  l)  PRINT  <5Q0*<NTITlFm,K«l.l3l 
r&TR7  uoo»i.J.<fM(.j.k1  TRm 7? J 
300  CONTINUE 

pfljMj  flofJ  —  —  -  — — - - - 

700  CONTINUE 

I^InEnD  ,LT.  10)  ltd  TO  t 
C  NFND  .GE.  10  WILL  HALT  PROGRAF 

Too — FwrrrtriTT - - — 

200  FORF AT ( a  1 1  o  I 

T075 - FTR-WiT(  A2.3rET2E17TrO) - 

300  FORF AT ( i H 1 1 

77To  EORNlTI  tNt.//i3A6. 

1  /3*. H-J.3X.lFK.6X.4F0000.12*.4F040Q»12X,4HOfl00.12X. 

r  •arreod/i ?y. «Hi'6T>'o«'igv. . gRtrrriT. bwpur? t - 

ROO  FORFAT(lHA) 

TTo73 — r0RVJTI2T4.6I2V.lPEi4.7U - 

ENC 


Flgura  Dl  Listing  of  PRIRTUTN,  an  algorithm  for  reading  and' 
printing  ths  1  MBs  atmospheric  noise  coefficients 


